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Abstract 

The neutrinoless double /3-decay of nuclei is reviewed. We dis- 
cuss neutrino mixing and 3x3 PMNS neutrino mixing matrix. Basic 
theory of neutrinoless double /3-decay is presented in some details. Re- 
sults of different calculations of nuclear matrix element are discussed. 
Experimental situation is considered. The Appendix is dedicated to 
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1 Introduction 



Observation of neutrino oscillations in atmospheric, solar, reactor and accel- 
erator neutrino experiments is one of the most important recent discovery in 
the particle physics. Small neutrino masses can not be naturally explained by 
the Standard Higgs mechanism. A new, beyond the Standard Model mech- 
anism of the generation of neutrino masses is required. The most plausible 
seesaw mechanism of the neutrino mass generation is based on the assump- 
tion of the violation of the total lepton number at a large scale and Majorana 
nature of neutrinos with definite masses. 

After it was established that neutrino masses are different from zero, 
the problem of the nature of neutrinos with definite masses z/j (Dirac or 
Majorana?) is the most actual one. Investigation of the neutrino oscillations 
can not allow to answer this fundamental question. The observation of the 
neutrinoless double /3-decay (Oi//3/3-decay) of some even-even nuclei would be 
a proof that z/j are Majorana particles. 

The neutrinoless double /3-decay is extremely rare process. First, this 
is a process of the second order of the perturbation theory in the Fermi 
constant. And, second, this process is possible due to helicity-flip. Thus, 
the matrix element of the process is proportional to the effective Majorana 
mass mi3j3 = U^iirii {nii is the mass of the neutrino z/j). Smallness of the 
neutrino masses is an additional reason for smallness of the probability of 
the Oz//3/3-decay. 

Very high values for the lower bounds of the half-lives of the Oz//3/3-decay 
of different nuclei were reached in the Heidelberg-Moscow [1], IGEX [2], 
CUORICINO [3] and other experiments. However, in order to reach the 
values of the half-lives of the Oz//3/3-decay which are expected on the basis 
of the neutrino oscillation data and the neutrino mass spectrum follows the 
inverted hierarchy, a new challenging experiments with a sensitivity to \mpp\ 
about two orders of magnitude better than the today's sensitivity are re- 
quired. It is expected that such a sensitivity will be reached in several future 
experiments. 

In this review we will consider in some details neutrino mixing. Then we 
will discuss the standard (Type I) seesaw mechanism of the neutrino mass 
generation. In the next section we will consider general properties of the 
neutrino mixing matrix and obtain its standard parametrization. Then we 
will discuss briefly the present status of neutrino oscillations. In the next 
section we will present quite detailed derivation of the matrix element of the 
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Oi//3/3-decay. Then we will consider effective Majorana mass under different 
assumptions about neutrino mass spectrum. In the next two sections we 
will discuss the present-day situation with the calculations of nuclear matrix 
elements of the Oi//3/3-decay and experiments on the search for neutrinoless 
double /9-decay. In the Appendix we will present a short biography of E. Ma- 
jorana and briefly discuss his 1937 paper in which the theory of the Majorana 
particles was developed and a possibility of the existence of such particles 
was discussed. 

For different aspects of the Oz//3/3-decay see reviews [H El EJ [TJ EJ HI [10] 

2 Neutrino mixing 

We will consider here the neutrinoless double /3-decay under two general 
assumptions 

1. The neutrino interaction is the Standard Model electroweak interaction. 
The Lagrangian of the standard charged current (CC) interaction has 
the form 

= --i^j^o^x)W-ix) + h.c. . (1) 

Here W^lx) is the field of the charged vector bosons, g is the 
constant of the electroweak interaction and 

Ja'^ix) = 2 Yl '^ILixhakix) + f^{x) (2) 
l=e,fi,T 

is the sum of the leptonic and hadronic charged current. The hadronic 
charged current is given by the expression 

= 2{uLixh^ drix) + Ci(x)7. srix) + hixha &r (x)), (3) 

where 

q=d,s,b q=d,s,b q=d,s,b 

(4) 

In dl]) the matrix V is the 3x3 Cabibbo-Kobayashi-Maskawa (CKM) 
quark mixing matrix [TTl [12] . 

The interaction ([1]) perfectly describe the data of numerous experiments 
on the study of the weak decays, neutrino reactions etc. 
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2. The neutrino mixing takes place. 

Neutrino fields viiix) in the leptonic current ([2]) are mixed fields 

3 



1=1 

Here i'i{x) is the field of neutrino with mass rrii and U is the 3x3 
Pontecorvo-Maki-Nakagawa-Sakata [T3l [H] neutrino mixing matrix. 

The hypothesis of the neutrino mixing was confirmed by the observation of 
the neutrino oscillations in experiments with the atmospheric, solar, reactor 
and accelerator neutrinos. All existing neutrino oscillation data are described 
if we assume that the number of massive neutrinos is equal to the established 
number of fiavor neutrinos (three). 

Quarks are charged particles; the quarks and antiquarks have the same 
masses and their charges differ in sign. Thus, the quark fields q{x) are com- 
plex Dirac fields. 

The electric charges of neutrinos are equal to zero. For neutrinos there 
are two fundamentally different possibilities. 

• If the total lepton number L = Lg + L^ + L^ is conserved, neutrino fields 
Vi{x) are complex four-component Dirac fields. In this case neutrinos z/j 
and antineutrinos Vi have the same mass and different lepton numbers 



• If there are no conserved lepton numbers, neutrino fields i'i{x) are two- 
component Majorana fields. In this case z/j = i/j. 

Investigation of the neutrino oscillations does not allow to distinguish these 
two possibilities [Ul ES]. In order to reveal the Majorana nature of i^j it is 
necessary to observe processes in which the total lepton number is violated. 
Neutrinoless double /3-decay of some nuclei is the only such process the study 
of which allows to reach the necessary sensitivity. 

The nature of neutrinos with definite masses and the form of the neutrino 
mixing is determined by the neutrino mass term of the Lagrangian. We will 
consider now possible mass terms for neutrinos (see [T7| IT8| [19j). 

A neutrino mass term is the Lorenz-invariant product of the left-handed 
and right-handed components of neutrino fields. The three left-handed cur- 
rent fields z//i(x), components of SU{2) doublets, must enter into any neu- 
trino mass term. If we assume that three right-handed singlet fields ^^{x) 




(5) 



= 1). 
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also enter into the Lagrangian in this case we can build the following neutrino 
mass term 

= -z/L(x)M°z/fi(x) + h.c., (6) 

where 



(7) 



and M° is the 3x3 neutrino mass matrix. It is obvious that the total 
Lagrangian with the neutrino mass term (E]) is invariant under the global 
gauge transformations 

Vl{x) e*^z/i(x), vr{x) e'^UR{x), Il,r{x) e'^lL,R{x), q{x) q{x), 

(8) 

where A is an arbitrary constant phase. The invariance under the transfor- 
mation dS]) means that the total lepton number L is conserved. 

The mass term can be easily diagonalized. For a complex matrix 
we have 

= U m V\ (9) 

where U and V are unitary 3x3 matrices and m is a diagonal 3x3 matrix 
{rriik = rrii 6ik, mj > 0). From and ([9]) we find 

3 

C^{x) = -u'^{x)mu'^{x) = -^m,z/i(a;)z/i(a;), (10) 



1=1 



where 




U^^L =\^2l\, Vr= V^Vl =\ ^2R \. (11) 




The expression (llUp is the sum of standard mass terms for the Dirac fields 
Vi{x) with masses m,. From ffTTl) we find that the flavor fields vwix) are 
connected with the left-handed components of the Dirac neutrino fields vnix) 
by the mixing relation 

3 

= J2UiiUiLix). (12) 
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We assumed that not only left-handed fields uil^x) but also the right-handed 
fields vir{x) enter into the total Lagrangian. In the original Glashow, Wein- 
berg and Salam papers |20l |2T1 [22], in which the Standard Model was pro- 
posed, it was assumed that only vil^x) fields, components of the lepton SU (2) 
doublets, enter into the Lagrangian. In the seventies after the success of the 
theory of the two-component neutrino it was natural to make this simplest 
assumption. In such a Standard Model with a SU{2) Higgs doublet neutrinos 
are massless particles. We can, however, generalize the original SM and to 
build a model in which neutrino masses and neutrino mixing are generated 
by the spontaneous violation of the symmetry in the same way as masses 
and mixing of quarks and leptons. In such a model the neutrino mass term 
is the Dirac mass term ([H]). 

We know from experiment that neutrino masses are many orders of mag- 
nitude smaller than masses of quarks and leptons. For example, for the 
particles of the third family 

mt - 173 GeV, ~ 4.2 GeV, ~ 1.78 GeV, mg < 2-10"^ GeV. (13) 

In the framework of the SM there is no natural explanation of such big differ- 
ence between masses of neutrinos and other fundamental fermions belonging 
to the same family. It is very implausible that small neutrino masses are 
generated by the SM Higgs mechanism. 

The small Dirac neutrino masses can be generated, however, in some 
models beyond the SM, for example, in the model with large extra dimensions 
[23] . In such a model the Newton law at small distances r has the form 
F = where n is the number of the extra dimensions and M is 

a new scale (~(1-10) TeV). Dirac neutrino masses in the model with extra 
dimensions are given by the expression 

TJii ~ kiV p. 

Here v ~ 250 GeV is the electroweak scale and (3 = ^ ^ (10^^^ - 10"^^) 
is a suppression factor (Mp ~ 1.2 10^^ GeV is the Plank mass). 

We will build now a neutrino mass term assuming that fields and 
i'ir{x) enter into the Lagrangian and there are no conserved lepton numbers. 
Let us consider the conjugated fields 

{i^Lr = c{vLr, {i^Rr = c{vRr, (m) 
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where C is the matrix of the charge conjugation which satisfies the relations 

C£C-' = -7„, = -C. (15) 

It is easy to show that [lyiY ((^k)^) is the right-handed (left-handed) com- 
ponent of the conjugated field. 

In fact, for the left-handed and right-handed components we have 

75 J^L = -l^L, 75 l^R = l^R- (16) 

From these relations we find 

i^Ll5 = T^L, i^Rl5 = -T^R- (17) 

Now, taking into account that C'-fJC"^ = 75, we have 

75 {i^lY = {r^LY, 15{i^rY = -{^rY- (18) 

From f lTS]) we conclude that {i^lY ^^"^ {^rY the right-handed and left- 
handed components. 

The most general neutrino mass term, which can be build from the flavor 
left-handed fields viLix) and sterile fields z/;/j(x)|l] has the form 

= -\ Mf{uLY - M° VR - ^ M^M^z/^ + h.c, (19) 

where columns vl,r are given by ([7]) and and M° are nondiagonal 

complex 3x3 matrices. It is easy to show that are symmetrical matrix. 
In fact, taking into account Fermi-Dirac statistics of the fields z/i,i?, we have 

-v,^n Ml^C-vl^ = (M,^^)^C^^I.^ = (M,^^)^ C^l^. (20) 

From this relation we find 

M^,n = {Mr,Rf. (21) 

It is obvious that the first and the third terms of the expression (fT9!) are not 
invariant under the global gauge transformations ul^r — )■ e^^u^j^. Thus, in 
the case of the mass term (fT9l) the total lepton number L is not conserved. 

^ Neutrino fields, which do not enter into the Lagrangian of tlie standard elecroweak 
interaction, are called sterile. 
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The first and the third terms of the expression (1191) are called the left- 
handed and right-handed Majorana mass terms, respectively. The second 
term is the Dirac mass term. The mass term £D+m jg ^g^ally called the 
Dirac and Majorana neutrino mass term [2H [25] . 

We will show now that in the case of the mass term ([T9l) neutrinos with 
definite masses are Majorana particles. 

The mass term C^+^ can be presented in the following form 



D+M 



iniM°+^K)^ + h.c. 



Here 



and 



riL 



M 



Mf 



(M 



(22) 
(23) 

(24) 



is a symmetrical 6x6 matrix. 

A symmetrical matrix M can be presented in the form 

M = UmU^, 



(25) 



where U is an unitary matrix and m is a diagonal matrix with positive 
diagonal elements. 

From ([22]) and ([25]) we have 



D+M 



1 



—III m 

-- V m V 



1 ^ 

4e 



m,; Vi. 



(26) 



Here 



From (12^ and (ITT)) we conclude that 



/ ^1 \ 



(27) 
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• The field z/j(x) {i = 1, 2, ...6) is the field of neutrinos with mass mj. 

• The field z/i(x) satisfies the Majorana condition 

v,{x) = vl{x) = CvJ{x). (28) 

Taking into account the unitarity of the matrix f/, from f l27|l we find 

nL = Uv^. (29) 



From fl29|l we obtain the following mixing relations in the general Dirac and 
Majorana case 

6 6 

m = ^ Uii UiL, [virY = ^ Uji UiL, (30) 

1=1 i=l 

where U is an unitary 6x6 mixing matrix and z/j is the field of the Majorana 
neutrino with mass mj. 

Let us discuss the meaning of the Majorana condition fl28p . A non her- 
mitian field z/(a;) can be presented in the following general form 

K^) = / ^ {crip) n^ip) e-'^^ + dlip) u^X-p) e^-) d^p, (31) 

where Cr{p) is the operator of absorption of neutrino with momentum p and 
helicity r and dl{p)) is the operator of creation of antineutrino with momen- 
tum p and helicity r and u^'{—p) = C{u'^{p))^. If the field z/(a;) satisfies the 
Majorana condition ( 128|) we find 

Cr{p) = dr{p). (32) 

Thus, if v{x) is the Majorana field, the neutrinos and antineutrinos are iden- 
tical particles. In other words the Majorana field is the field of truly neutral 
particles. There is no notion of particles and antiparticles in the case of the 
Majorana fieldl 

We will finish this section with the following remarks 

1. Dirac and Majorana mass term can be generated only in theories be- 
yond the SM. 



^In the case of the Majorana field there are no conserved charges which ahow to dis- 
tinguish particles and antiparticles. 
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2. If we assume that only left-handed fields enter into the mass 
term and the lepton number is not conserved, we come the following 
(Majorana) mass term|26] 

C'^ = -\v-LMf{vLr + \..c., (33) 

where is 3 x 3 symmetrical matrix. After the diagonalization, the 
mass term (133!) takes the standard form 

1 ^ 

C"^ = --Y^^i^i^^ (34) 

1=1 

and we come to the Majorana mixing 

3 

i'iL = Y^Uiii'iL. (35) 

1=1 

Here [/ is 3 x 3 mixing matrix and Vi is the Majorana field with the 
mass nii which satisfies the condition ( 128|1 . Notice that Higgs triplet is 
needed for the generation of the mass term 

3. From the Majorana condition ( 128|) we have 

= = M'- (36) 

Thus, in the case of the Majorana field right-handed and left-handed 
components are connected by the relation (!36|) . In the case of the Dirac 
field right-handed and left-handed components are independent. This 
is the major difference between Majorana and Dirac fields. 

Right-handed components of neutrino fields enter into the Dirac mass 
term. If neutrinos are massless there are no mass term in the La- 
grangian. This is the reason for the well known theorem [27] which 
states that it is impossible to distinguish massless Dirac and Majorana 
neutrinos in the case of left-handed interaction. 

4. The Dirac and Majorana mass term opens a possibility of the exis- 
tence of the sterile neutrinos. If masses are small, in this case in 
addition to the mixed flavor left-handed neutrinos u^, and Vr mixed 
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left-handed antineutrinos un, quanta of mixed right-handed fields ui^, 
must exist. Because right-handed fields do not enter into the standard 
CC and NC interactions, have no electroweak interaction. They 
are called sterile neutrinos. Let us notice that the existed LSND indi- 
cation in favor of the sterile neutrinos [28] was not confirmed by the 
MiniBooNE experiment |29]. 

5. In the case of the Dirac and Majorana mass term there are additional 
sterile right-handed fields uui and many parameters in the mass ma- 
trix. This mass term open a possibility to explain the smallness of the 
neutrino masses. This (so called seesaw) possibility will be considered 
in the next section. 

3 Seesaw mechanism of the neutrino mass 
generation 

The most popular mechanism of the generation of small neutrino masses 
is the seesaw mechanism [30]. In order to explain the main idea of this 
mechanism we consider the simplest case of one generation. The Dirac and 
Majorana mass term is given in this case by the expression 



We will assume that mL^mc and mn are real parameters. Let us write Eq. 
fl37|) in the matrix form. We have 





(37) 




(38) 



Here 




(39) 



and 




(40) 



The real symmetrical matrix M 



can be presented in the form 
= 0m'0^, 



(41) 
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where 

— sin 6 cos d 

and m^^ = m^6ik, m'^ being an eigenvalue of the matrix M^~^^. We have 



0= "1. 1 (42) 



I , . 1 



m 



1.2 



- {ttir + mi) =F 2 V i'^R - + 4 m|,. (43) 



From (HT]) . P2|) and fH3|) for the mixing angle 6' we obtain the following 
relations 

tan2t^= , cos2t^ = — = . (44) 

niR - niL \/{mR - miY + Arri^) 

The eigenvalues m'^ ^ can be positive or negative. Let us write down 

= mir]i , (45) 

where = |m'-| and rji = ±1. 
From dH]) and (jlS]) we find 

M^^"^ = UmU^ , (46) 

where 

f/ = 0y7/ (47) 

is an unitary matrix. Using the general results of the previous section, we 
easily bring the mass term (l48ll to the standard form 

i=l,2 

Here 

= U^UL + (f/^nz.)^ = ( ) ' ^^^^ 
z/j being the Major ana field with the mass mi. From (H9|) we have 

riL = f/ . (50) 
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Thus, the fields pl and {ur)^ are connected with the fields un and z/2l by 
the following mixing relations 



vl = cos 6y/r]l viL + sin 9.Jr]i V2L 
{vrT = -sin6'v^z/iL + cos6'yr7^z/2L (51) 

Neutrino masses are many orders of magnitude smaller than masses of lep- 
tons and quarks which are generated by the standard Higgs mechanism of 
the electroweak symmetry breaking. This fact is commonly considered as 
an evidence in favor of a non-standard mechanism of neutrino mass genera- 
tion. The seesaw mechanism connects smallness of neutrino masses with the 
violation of the total lepton number at a very large scale. 

The standard (type I) seesaw mechanism [30] is based on the following 
assumptions 

1. There is no left-handed Majorana mass term in the Lagrangian [rriL = 
0). 

2. The Dirac mass term is generated by the Higgs mechanism {m£, is of 
the order of the mass of a charged lepton or quark). 

3. The constant m^, which characterize the right-handed Majorana mass 
term, the source of the violation of the total lepton number, is much 
larger than rriD: 

niR > iriD (52) 

From (is]), ^ and (ES]) we have 

mi ~ mo <^ — mji, tanu ~ ^ 1. (53) 

rriR tur 

Thus, the seesaw mechanism generates Majorana neutrino mass mi which 
is much smaller than a Dirac mass of a lepton or quark. As a consequence 
of the seesaw mechanism a heavy Majorana particle with a mass m2 — mR 
must exist. 

Let us consider now the case of the three families. The seesaw mixing 
matrix has in this case the form 
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Here m/) is a complex 3x3 matrix, is a symmetrical complex matrix 
and rriD Mr. 

Let us introduce the matrix M by the relation 

jjT ^seesaw JJ^^, (55) 

where U is an unitary matrix. We will show now the matrix U can be chosen 
in such a form that M is the block- diagonal matrix. 

Notice that in the case of one generation up to terms linear in — ^ 1 
we have 

U^'^ - (n^, 7 ■ (56) 



Let us consider the matrix 



C-(_^, f). (57) 

where A is a 3 x 3 matrix and Aik ^1. It is easy to see that up to linear 
in A terms WU 1. The non-diagonal element of the symmetrical matrix 

over A approximation is equal to 

Mr Al (58) 



If we choose 

A^ = M^^ ml (59) 
the matrix f/^ ]\^seesaw jj ^g^j^gg block- diagonal form 

jjT ^seesaw f/ ^ |^ "^^^ ^ . (gQ) 

For the left-handed Majorana neutrino mass term from (160!) we find 

C^'^ = -\ yLMf{ur.r + \i.c. (61) 

where 

Mf = -rriD M^^ ml (62) 

and is given by ([7]). 

The Eq. fl6T|) is the mass term for three light Majorana neutrinos. After 
the diagonalization of the total mass term in addition to Majorana neutrino 
mass term we will obtain a mass term for three heavy Majorana particles. 
Thus, in the case of the Dirac and Majorana mass term with the matrix ([54 
in the spectrum of masses there are 
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• Three light Majorana neutrino masses. 

• Three heavy Majorana masses, which are characterized by the scale of 
the violation of the total lepton number. 

These are general features of the seesaw mechanism. The values of neutrino 
masses and mixing angles can be obtained only in the framework of a concrete 
model. 

Thus, the seesaw mechanism connects smallness of the neutrino masses 
with violation of the total lepton number at a large scale|f| The observation 
of the neutrinoless double /3-decay would be an evidence in favor of this 
mechanism. Let us notice that the existence of heavy Majorana particles, 
seesaw partners of neutrinos, could allow to explain the baryon asymmetry 
of the Universe (see [31] ) . 



4 Neutrino mixing matrix 

In this section we will consider the general properties of the unitary 3x3 
Dirac (or Majorana) mixing matrix. 

An unitary n x n matrix U is characterized by ra^ real parameters]^ The 
number of the angles which characterize the unitary n x n matrix coincides 
with the number of parameters which characterize a real orthogonal n x n 
matrix O {0^0 = 1). Thus, for the number of the angles we hav^ 

^(^-1) 

^ang = . (63) 

Other parameters of the matrix U are phases. The number of phases is equal 
to 

9 n(n — 1) n(n + 1) , 
^Ph = n'- = (64) 

The number of physical phases in the neutrino mixing matrix is smaller than 
nph. The neutrino mixing matrix enters into the charged current. Let us 



■^Usually it is assumed that this scale is about (10^^ — 10^^) GeV. 

^In fact, it can be presented in the form U — e*^, where H is the hermitian matrix. 
The hermitian matrix is characterized by n + 2 ( " — -n? real parameters. 

^The orthogonal matrix O can be presented in the form O = e^, where ~ —A . 
Diagonal elements of the matrix A are equal to zero. The number of the real non-diagonal 
elements is equal to . 
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consider first the case the of Dirac neutrinos z/j. Because phases of the Dirac 
fields and vnix) are arbitrary, the matrices U and 

U' = S\f3) U S{a) (65) 

are equivalent. Here Sii'{/3) = e^^^dui, Sik{oi) = e^°'^6ik and (3i, at are real, 
arbitrary phases. 

We can use this freedom in order to exclude {2n — 1) phases from the 
matrix f/J§ Thus, in the case of the Dirac neutrinos the number of the 
physical phases in the mixing matrix U is equal to 

n(n + l) , , (n-l)(n-2) , , 

^Ph = - (2n - 1) = ^ \ ^. (66) 

In the case of the mixing of the three Dirac neutrinos the mixing matrix is 
characterized by three mixing angles and one phase. 

Let us consider now the case of the Majorana neutrinos Vi. The Majorana 
condition 

vlix) = u,ix) (67) 

does not allow to include arbitrary phases into the Majorana fields. For the 
number of the physical phases we have in the Majorana case [151 [IS] 

nX = '-^-n='^^. (68) 

Thus, in the case of the three Majorana neutrinos the mixing matrix is char- 
acterized by three mixing angles and three phases. 

We will obtain now constraints on the neutrino mixing matrix which 
follow from the condition of the CP invariance in the lepton sector. Let us 
consider first the Dirac neutrinos z/j. The condition of the CP invariance in 
the lepton sector has the form 

VcP £r(x) = £f^(x'). (69) 

Here Vcp is the operator of the CP conjugation, x' = (x°, —x) and 

(70) 



^ We can always make one element of the matrix S{a) (or S{f3)) equal to one. In fact, let 
us present the matrix S{a) in the form S{a) — e'"" S{a), where ai — ai — an- The phase 
factor e'"" can be, obviously, included into S'^(/3). We have in this case e*"" — 5^(/3), 

where Pi = Pi - an- 
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is the Lagrangian of the CC interaction of neutrinos, leptons and W^-bosons. 
Taking into account arbitrariness of the phases of fermion fields, we can put 
CP phase factors of the lepton and neutrino fields equal to one. We have 

VcplL{x)Vc^ = ^'Cmx'), Vcpv,l{x)VcI. = i' C vUx'). (71) 

From these relation we find 

VcP Il{x) = -ll{x') VcP v^L{x) Vcl = -^Ux) C-V- (72) 

For the field of the charged vector bosons we have 

VcP W^{x) Vcl = -6a Wlix'), (73) 

where 5a is a sign factor {6q = 1, 6i = —1). From all these relations we easily 
find 

VcpCf'^{x)V^^ = -^J2^^L{x')^aUulL{x')W''{x') 

- j^J^lL{x')iaU;,u,,{x')W'^Hx'). (74) 

Comparing (169!) and fl7^ we come to the conclusion that in the case of the 
CP invariance in the lepton sector the Dirac mixing matrix is real 

Uh = U;,. (75) 

We will consider now the case of the Majorana fields [521 ESI El]- The CP 
transformation of the Majorana field z/j has the form 

VcP Ui{x) Vch = v: 7° C vj{x') = T]* 7° ^^(x'), (76) 

where t]* is a phase factor. Unlike the Dirac fields, it can not be included in 
the field. We will show now that the phase factor rji can take the values ±i. 
In fact, from fl76l) by the hermitian conjugation and multiplication from the 
right by the matrix 7*^ we find 

Vcpu,{x)Vc^ = V^Mx')l''■ (77) 
From this relation we have 

VcP C yj{x) = mC 7°^ C uf{x') = -V, 7° C uf{x'). (78) 



17 



Finally, taking into account the Majorana condition we find 

VcP v.^x) = 7° (79) 
If we compare now f l76l) and fl79|) we conclude that 

V: = -Vi. ^' = -1- (80) 

Thus, the CP parity of a Majorana field can take values ±z. 

From flB^ . f lTT]) and fl7^ we find that in the case of the CP invariance in 
the lepton sector the Majorana mixing matrix satisfies the condition 

Uu V* = u;,. (81) 

Finally we will obtain the standard parametrization of the 3x3 Dirac 
mixing matrix. Let us consider two systems of orthogonal and normalized 
vectors \i) and l^i) {i = 1,2,3, / = e,fi,T). We have 

{k\t) = S,k, m = Si'i. (82) 
Vectors {ui) and \i) are connected by the relation 



From ( !82|) it is obvious that U is an unitary matrix. 

In the most general case vectors |z/;) can be obtained from vectors \i) by 
three Euler rotations. The first rotation will be performed at the angle 612 
around the vector |3). New orthogonal and normalized vectors are 

=Ci2 \l)+su |2) 

|2)(1) =-Si2 |1)+C12 |2) 

|3)W =|3), (84) 



where cu = cos ^12 and S12 = sin6'i2. In the matrix form (|8^ can be written 
as follows 

|zy)(l) = f/(l) \iy). (85) 

Here 

= I |2)(i) I , Iz.) = I |2) I (86) 
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and 



Cl2 Si2 
-Si2 Ci2 

1 



(87) 



We will perform now the second rotation at the angle 613 around the vector 
1 2)*^^). At this step we will introduce the CP phase 5. We have 



Here 




(89) 



Finally, let us perform the rotation around the vector ll)^^-* at the angle ^23- 
We have 

l^mix^ = f/{3) |^){2)^ (gQ) 



Here 




(91) 




and 

From dHS]), dHH]) and dlO]) we find 

Iz/'"") =U* 

where 



(92) 



U = (t/(3) f/(2) jjiDy 



1 

C23 S23 
— S23 C23 




Si3e 




(93) 



C12 S12 
-S12 C12 1(94) 
1 
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This is so-called the standard parametrization of the 3x3 Dirac mixing ma- 
trix. This matrix is characterized by three mixing angles 612, O23 and 613 and 
the CP phase 5. From flMl) we have 

(C13C12 C13S12 si3e""^ \ 

-C23S12 - S23Cl2Sl3e"' C23C12 - S23Si2Si3e*'' C13S23 • (95) 

iS iS I 

-523512 ~ C23Ci2-Si3e — -S23C12 ~ C23Si2-Si3e C13C23. / 

The 3x3 Majorana mixing matrix is characterized by three mixing angles 
and three CP phases. It can be presented in the form 

^ jj gM^^^^ (gg) 

where the matrix U is given by flMl) and 



(a) = I e^°2 I ^ (97) 
where ai,2 are additional Majorana phases. 



5 On neutrino oscillations 

The most important manifestation of the neutrino mixing are neutrino os- 
cillations. Neutrino oscillations are based on the fact that in processes of 
neutrino production and neutrino detection due to Heisenberg uncertainty 
principle small neutrino mass-squared differences can not be resolved. As a 
result, in a weak decay 

a^b + l^ + ui (98) 

together with the lepton a "mixed" left-handed flavor neutrino ui is pro- 
duced. The state of ui is a coherent superposition of the states of neutrinos 
with definite masses 

h) = Y.^uW^), (99) 

i 

where |z/j) is the state of neutrino with mass rrii and momentum pi. 

If at t=0 flavor neutrino z/; is produced, at the time t for the neutrino 
state we have 

3 3 
kO* = e-^^*|z.,) e-^^-*f/* h)=J2\^i')J2Ui'^e-'''^'U;,. (100) 

i=l I' i=l 
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Thus, the probabihty of the transition ui — )■ z/;/ during the time interval t is 
given by the expression 



P{ui -> u,) = I J]t/,ie-^(^-^'=)*f/;J^ (101) 

i=l 

where k is fixed. If all phase differences are small {\Ei — Ek\t 1) or/and 
there is no mixing {U^ = du) in this case it will be no neutrino oscillations 
{Pi^i — )■ vv) — Sill). Thus, neutrino oscillations are effect of the neutrino 
mixing and relatively large phase difference(s). 

Assuming that pi = p, we obtain the standard expression for the phase 
difference 

[E, -E,)t^^L. (102) 

Here Am^^ = mf — m\ and L ~ t is the distance between a neutrino source 
and neutrino detector. 

The transition probability P{vi — t- vi') depends on six parameters (two 
mass-squared differences Amga and Am^2) three mixing angles ^23, ^12 and 61^ 
and CP phase 5). However, from analysis of the data of neutrino oscillation 
experiments follows that the parameter sin^ ^13 and the ratio are small: 

^^3■10-^ sin2^^i3<5-10-^. (103) 

If we neglect contribution of the small parameters to the transition prob- 
abilities, we will find that in the atmospheric and accelerator long baseline 
region of the values of the parameter ^ the two-neutrino ^ Ur oscillations, 
driven by Am|3, take place. From f llOip and f ll02p for the probability of 
to survive we obtain the following expression (see, |18] ) 

P(z/^ ^ z,^) ~ 1 _ i sin^ 2^23 (1 - cos Am^g A) . (104) 

In the KamLAND region v^. ^ v^^j. oscillations, driven by Am^g! ^ake place 
in the leading approximation. For the probability of z/g to survive we obtain 
the following expression (see [T8] ) 

V{ve -^Ve)^l~]: sin^ 2^^i2 (1 - cos Aml^^) . (105) 
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In the leading approximation the probabihty of solar u^, to survive in matter 
is also given by the two-neutrino expression. It depends on tan^6'i2, ^rnl^ 
and electron number density in the sun. 

We will present now the results of the analysis of the experimental data. 
From the analysis of the data of the atmospheric Super- Kamiokande exper- 
iment for the parameters Am23 and sin^ 26*23 the following 90 % CL ranges 
were obtained [53] 

1.5 ■ 10"^ < Am^3 < 3.4 ■ 10" sin^ 2^23 > 0.92. (106) 

The results of the atmospheric Super-Kamiokande experiment were con- 
firmed by the K2K |36j and MINOS accelerator long-baseline neutrino 
oscillations experiments. From the analysis of the MINOS data for the neu- 
trino oscillation parameters the following values were found [37] 

Am^3 = (2.43 ± 0.13) ■ 10"W, sin^ 2^23 > 0.90 (90% CL). (107) 

From the global analysis of the data of the reactor KamLAND experiment 
and data of the solar neutrino experiments for the parameters Am^2 
tan^ 6i2 the following values were obtained [55] 

Aml^ = ■ 10-^ eV^ tan^ 6^ = OA7tom (108) 

In the reactor CHOOZ experiment no indications in favor of — ?• Ug 
transitions, driven by Am23, were found. From the exclusion plot, obtained 
from the data of this experiment, for the parameter sin^ ^13 the following 
upper bound can be inferred 

sin2ei3<5-10-2 . (109) 

At present a stage of the high precision neutrino oscillation experiments 
starts. In the future DOUBLE CHOOZ [ID], Daya Bay p[T] and RENO [12] 
reactor neutrino experiments sensitivities to the parameter sin^ 26*13 will be 
(10-20) times better than in the CHOOZ experiment. The same sensitivity is 
planned to be reached in the accelerator T2K experiment [13|. In this exper- 
iment parameters Am23 ^^^^ ^6^23 will be measured with the accuracies 
^(Amlg) ~ 10~^eV^ and 5(sin^ 26^23) ~ 10~^, correspondingly. High precision 
neutrino oscillation experiments are planned at the future Super Beam [H], 
Beta-beam [4S], and Neutrino Factory facilities [15] . 
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6 Basic elements of the theory of Oz//3/3-decay 

In this section we will consider the neutrinoless double /3-decay of even-even 
nuclei [H |5] 

{A,Z)^{A,Z + 2) + e- + e-. (110) 

We will assume that 

• The Hamiltonian of the weak interaction is given by the SM. 

• The neutrino mixing takes place. 

• Neutrinos with definite masses Ui are Majorana particles. 
For the effective Hamiltonian of the process we have 

ni{x) = ^2 ^ CLixh^ Uu PiL{x) rix) + h.c. (Ill) 

* i 

Here is the Fermi constant, j"(a;) is the hadronic charged current and 
the field Uii^x) satisfies the condition 

u^{x) = Cufix) = iy,{x). (112) 

The neutrinoless double /3-decay is the second order in Gp process with the 
virtual neutrinos. The matrix element of the process is given by the following 
expression 

(/l^'K) =4^!!! [^Jn,,N,,I Y^uMe''^'"'!. f/e.(0|r(z.,i(xO uf^{x,)\0) 

X7j Ue^ul{p2)e'P'''^Nf\T{J"{x,)J^{x2))m d^xid'x2 - {pi ^ P2) ■ (113) 

Here pi and p2 are electron momenta, J"(x) is the weak charged current in 
the Heisenberg representation, Ni and Nj are the states of the initial and 
the final nuclei with 4-momenta Pi = {Ei,pi) and Pj = {Ef,p'f), respectively, 
and N„ = ^ , — is the standard normalization factor. 

^ (27r)3/2y^ 

Let us consider the neutrino propagator. From the Majorana condition 
f im we find 

{0\T{u,L{x,)uJ^ix2m = -^-^{0\T{u,{x,Mx2m C. (114) 
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Further, we have 

(0|T(z.,(^0^.(^2))|0) = 7^ / e-^'^ (115) 

Thus, we for the neutrino propagator we find the following expression 

{0\TMx^)u,Lix,)m = J e-^^ ^^^"^^^^^ A C. 

(116) 

The neutrino propagator is proportional to rrii. It is obvious from flllSp that 
this is connected with the fact that only left-handed neutrino fields enter into 
the Hamiltonian of the weak interaction. In the case of massless neutrinos 
{rrii = 0, = 1,2,3), in accordance with the theorem on the equivalence of 
the theories with massless Majorana and Dirac neutrinos, the matrix element 
of the neutrinoless double /3-decay is equal to zero. 

Let us consider the second term of the matrix element (I113p . It is easy 
to show that 

UliPlhai^ - l5hl3Cul{p2) = Ml(P2)C^7J(1 " llh^uKpi) 

= -UL{p2hl3{l - IbhaCulipl) . (117) 

If we take into account flll7p and the relation 

T{J^{x2).r{x,)) = T{r{x,)j\x2)) (118) 

we can show that the second term of the matrix element flllSp is equal to the 
first one. Thus, for the matrix element we obtain the following expression 

(/l^'N) = -\-^) A^P.Ar,,y^,(pOe-^^--7.^E^>-y ^2 _ A 

x^^7/3C ul{p2)e'^'''HNf\T{r{x,)j\x2))\m d'^x^d^x^ (119) 

Initial nuclei in the process f lllOp are ^^Ge, ^^^Xe, ^^°Te, ^°°Mo and other 
heavy nuclei. The calculation of the nuclear part of the matrix element of the 

^Notice that in the case of the Dirac neutrinos {{)\viL{xi)vfj^{x2)\Q) = 

{Q\i'i{xi)vJ {x2)\Q)^—^ — 0. The neutrinoless double /3-decay is obviously forbidden 
in the Dirac case. 
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Oi//3/3-decay is a complicated nuclear problem. In such a calculation different 
approximations are used. We will present now the matrix element of the 
Oi^/3/3-decay in a form which is appropriate for such approximate calculations. 

Let us perform in flllQp the integration over the time variables and x^. 
The integral over X2 can be presented in the form 

...dxl = / ...dxl+ / ...dxl . (120) 

-00 J —00 J x1 

After the integration over in the neutrino propagator, in the region x^ > 
we finely 

iq {X1-X2) 1 r p-iQi (x°-x'^)+iq (X1-X2) 

A = TTTT^ / TT-a d'q , (121) 



(27r)4 J q^-mf ^ (27r)3 j 2 q] 

where 



In the region x^ < x\ we have 

7-^ / —o = T-^ / n d^Q- (123) 

(27r)4 J q^-mf ^ (27r)3 J 2 q^ ^ I ^ 

For the operators J"(x) from the invariance under the translations we have 

r{x) = e'^"V"(x)e-'^"°, (124) 

where H is the total Hamiltonian. From this relation we find 

{Nf\rix^)J^ix2m) = 

J2e^i^f-^^Ke^i^--^^H (iV;|J-(fl)|A'n)(iV„|J^(x2))|iV,) , (125) 

n 

where |A^„) is the vector of the state of the intermediate nucleus with 4- 
momentum P„ = {En,]fn). In f ll25p the sum over the total system of the 
states |A^„) is assumed. 

Taking into account that at ±oo the interaction is turned off we have 

/ 6^"^^° dxl / e^("-*^)^2 dxl = lim (126) 

J-oo J-oo a-te 



oo 

is assumed that in the propagator mf = m| — ie 
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and 

e^'^^°rf4^/ e*('*+^^)^° rfx^ = lim . (127) 

^0 Jo ' e^oa + ie 

From ffT26D and ffT27l) we find 



I 



,.,?<4-.!) . V (iV/|J°(a)|iV„>(A-.|J^(f.))|iV, + po ^ po _ j;.) (123) 

^ -C'n ~^ P2 ~^ Ei 16 

Taking into account all these relations, for the matrix element of the neutri- 
noless double /3-decay we obtain the following expression 



1 f e*^^(-^i--^) ,3 (iV^|J°(a:-l)|iV„)(iV„|J'5(xl))|iV,) 



^ ^ (iv.|.^^.)|iv )(iv„|.-(.-0)|iv.> 

^ En+Pi + ql - Ei-'>^ 



The equation (11291) is the exact expression for the matrix element of Ou^/S- 
decay in the second order of the perturbation theory. We will consider major 
0+ —7- 0+ transitions of even-even nuclei. For such transitions the following 
approximations are standard. 

1. Small neutrino masses can be safely neglected in 

The averaged momentum of the virtual neutrino is given by the relation 
g ~ i, where r is the average distance between two nucleons in nucleus. 
Taking into account that r ^ 10~^^ cm, we have q ~ 100 MeV. Neutrino 
masses are smaller than 2.2 eV. Thus, we have q^ = + mf ~ q 

2. Long- wave approximation. 

We have PkXk < PkR^ where R ~ 1.2 A^/^ ■ 10"^^ cm is the radius of 
nucleus {k = 1,2). Taking into account that Pk ^ ^ Mev, we have 
PkXk ^ 1- Thus, e~*^^^^~*^2^'2 ^ 1 i.e. two electrons are emitted in 
5'-states. 
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3. Closure approximation. 

Energy of the virtual neutrino is much larger than the excitation energy 
{En — Ei). Thus, we can change the energy of the intermediate states 
En by average energy E. In this (closure) approximation we have 

{Nj\J-ix,)\Nn){Nn\J^ix,))m ^ {Nf\J-ix,)J'^ix-'2))m 

En+vl + q^i - Ei-ie E + pl + q - Ei - le ' 

4. The impulse approximation for the hadronic charged current J°(x). 

Taking into account the major terms, the hadronic charged current 
takes the form |^ 

r{x) ^ ^ 5(x - rn) rl[ gv{q^)g'^' + ^/^(g'Xr] • (131) 

n 

Here gv{q^) and gA{q^) are vector and axial formfactors , cTj and Xj are 
Pauli matrices, t+ = | (ti + 2x2) and index n runs over all nucleons in 
a nucleus. We have gviS^) = 1, 5'a(0) = gA — 1-27. 

It is obvious that r" = 0. Thus, in the impulse approximation the hadronic 
currents satisfy the relation 

J"(Xi) J^{X2) = J^{X2) J"(fi). (132) 

Further, the matrix 70,7^ in the leptonic part of the matrix element (I129p 
can be presented in the form 

7a7/3 = gal3 + ^ (7a7/3 " 7/37a)- (133) 

It follows from fll32p that the second term of (11331) does not give contribution 
to the matrix element. From (I13ip we have 

J'^ix-DUx-h) = 5^r^r-5(xi - Vn) 5{x2 - f^KgUQ^) ' dW) " , 

(134) 

^The pseudoscalar term in the one-nucleon matrix element of the hadronic charged 
current induces a tensor term in the current. From numerical calculations follow that its 
contribution to the matrix element can be significant (see 1471). 
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Neglecting nuclei recoil, we obtain in the laboratory frame 

where Mj and Mf are masses of the initial and final nuclei. From this relation 
we find 

,^A,^E-M,^,±(^l^)^E^Ml±Ml (135) 

The term ( ^^^^^ ), is much smaller than all other terms in the right-hand side 
of this relation. Neglecting this term, we have 

q + p° 2 + E - M,c^q + E (136) 



Further, taking into account that gy[<f) — 2 and gAiq 



0.71 GeV2 

a2 ^2 



where M4 1 GeV , we can neglect g -dependence of the formfactors. After 
the integration in the matrix element (I129p over xl and X2, for the neutrino 
propagator we find the following expression 

^ Hir^m.E), (137) 



(27r)3y g(g + ^_i(Mi + M/)) 47ri? 
where 

H{r,E) = ^-^ r-_singrrfg _ 

7rr Jo g + ^-|(M, + M/) 

Here i? is the nuclei radius and fnm. = — Vm- 

Taking into account all these relations, from fll29p for the matrix element 
of Oz//3/3-decay we obtain the following expression 

xM"" 5{p\ + pl + Mf-Mi), (139) 

where 

mpp = Y,Ulmi (140) 

i 

is the effective Majorana mass and 

9a 
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is the nuclear matrix element. Here 



= {^f\J2H{r,,,^,E) T^r^\^^) (142) 



is the Fermi matrix element and 

M'gt = {^f\Y,HK,^,E) tIt"^ a" ■ a"^)|^,) (143) 

n,m 

is the Gamov- Teller matrix element. In ( I142p and (11431) are wave 

function of the initial and final nuclei. 

From (I139P we conclude that matrix element of Oz//3/3-decay is a product 
of the effective Majorana mass m^^, the electron matrix element and the 
nuclear matrix element which includes neutrino propagator (neutrino poten- 
tial). Taking into account that E — \{Mi + Mj is much smaller than q for 
the neutrino propagator we obtain the following approximate relation 

, 2R p singr^ R 

Hir) ~ — / —dq = — . (144) 

n Jq qr r 

Using the standard rules, from (11391) we can easily obtain the decay rate of 
the Oi//3/3-decay. The electron part of the decay probability is given by the 
trace 

Tr(l + 75) (7 ■ p2 - mc)(l - 75) (7 ■ pi + me) = 8piPi . (145) 

Taking into account the final state electromagnetic interaction of the elec- 
trons and nucleus for the decay rate of the Oi//3/3-decay we find the following 
expression 

dV'^ = |m^^|2|M°-p^G^^-1^7l(EiE2-piP2Cose)x 

F{Ei, {Z + 2)) F{E2, {Z + 2)) p.p^ smOdO dE^, (146) 

where Ei^2 = Pi,2 is electron total energy {E2 = Mi — Mj — Ei), 9 is the angle 
between electron momenta pi and p2 and 

is the Fermi function {1] = Za ^). 
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From (11461) follows that for the ultra relativistic electrons ^-dependence 
of the decay rate is given by the factor (1 — cos^). Thus, ultra relativistic 
electrons can not be emitted in the same direction. This is connected with 
the fact that the helicity of the high energy electrons, produced in the weak 
interaction, is equal to -1 . If electrons are emitted in the same direction, the 
projection of their total angular momentum on the direction of the momen- 
tum is equal to -1. It is obvious that such electrons can not be produced in 

— )■ O'^ transition. 

From expression fll46p for the total decay rate we obtain the following 
expression 

= ;i = M \M'^? G^^iQ, Z), (148) 

^1/2 

wher^ 

G'^'^iQ, Z) = G% g\ £ dT^ £ sin OdO {E.E^ - pm cos 9)p,p, x 

F{E,,{Z + 2)) F{E,,{Z + 2)). (149) 

Here Ti = Ei — rrie, Q = Mi — Mf — 2me is the total released kinetic energy 
and is the half-life of the In the Table I we present numerical 

values of G^'^{Q,Z) for some nuclei [i8] . 

Table I 



The values of the factor G^'^{Q, Z) for some nuclei 



Nucleus 


G'^''{Q,Z) in units lQ-''^y-^eV-'' 


76Ge 


0.30 


looMo 


2.19 


130 rjg 


2.12 


i36Xe 


2.26 



The total rate of the Oz//3/3-decay is the product of three factors: 

1. The modulus squared of the effective Majorana mass. 

2. Square of nuclear matrix element. 

""^"An additional factor i is due to the fact that in the final state we have two identical 
electrons. 
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3. The known factor G°^(g, Z). 

We have considered in some details neutrinoless double /3-decay of nuclei 

{A,Z) {A,Z + 2) + e- +6-^. (150) 

There could be other second order in the Fermi constant Gp processes with 
the virtual Majorana neutrinos in which the total lepton number is changed 
by two. The examples are the decays 

K-^n+ + lx-+e- (151) 

and 

^ TT" + /i+ + /i+ (152) 

the process 

fx- + {A,Z) ^ {A,Z -2) + e+ (153) 

and others. 

The leptonic part of the operator which give contribution to matrix ele- 
ments of f ll50p . (115 ip and other similar processes is given by 

5^r(ri(xi)7a Uu (0| T{UiL{Xi)ulLix2)) \0) U,i -f^ T[ (x^)) , 1,1' = 6, fi, 

(154) 

where Majorana neutrino propagator is given by the expression flll6p . Taking 
into account that mf <^ p^, we can neglect in the denominator of the 
propagator. Thus, the matrix element of a process in which a lepton pair 
(//') is produced, is proportional to 



mil' 



^UiiUi'imi. (155) 



Analogously, matrix elements of the processes (I152p . (I153P and other similar 
processes are proportional to m^,. 

The sensitivities to the parameter \miii\ of the experiments on the search 
for the processes f ll5ip . fll52p . f ll53p and other similar processes are much 
worse than the sensitivity of the experiments on the search for Oz//3/3-decay 
to the parameter |m/3/3|. 

For example, in the experiment [52] on the search for the process fi^Ti 
e~^Ca. the following upper bound was obtained 

r(/i~Ti — )• e"'"Ca) in-12 



r(/i-Ti ^ all) 
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<1.7-10"'^ (156) 



For the probability of the decay — ^ vr ^~^fi~^ the following upper bound 
was reached [53]: 

T(K+ n-n+a+) „ , , 

From these data the following upper bounds can be found (see [7]) 

|m^e|<82MeV; |m^^| < 4 ■ 10^ MeV. (158) 

These values must be compared with the sensitivity of the experiments on the 
search for Oi//3/3-decay to the effective Majorana mass (in today's experiments 
\m0p\ ~ (0.2 - 1.3) eV (see below)). 

The effective Majorana mass is determined by neutrino masses and neu- 
trino mixing angles. An information about the neutrino mixing angles dik 
and neutrino mass-squared differences ^rnjf. was obtained from the data of 
the neutrino oscillation experiments. Taking into account these data, we will 
consider now possible values of the effective Majorana mass. 



7 Effective Majorana mass 

From neutrino oscillation data follows that one mass-squared difference (so- 
lar) is much smaller than the other one (atmospheric). For three massive 
neutrinos two types of neutrino mass spectra are possible in this case. 

1. Normal spectrum 

nil < 1712 < 'ms', Am^2 ^ ^^23 (159) 

2. Inverted spectrunj"] 

ms < nil < 1712] Am^2 ^ I^^^Ll (160) 

^^In order to have the same notation Amfj for the solar-KamLAND neutrino mass- 
squared difference and to determine this quantity as a positive one the neutrino masses are 
usually labeled differently in the cases of the normal and inverted neutrino mass spectra. 
In the case of the normal spectrum Am^^ > and in the case of the inverted spectrum 
Amf3 < 0. Thus, with such a notation the character of the neutrino mass spectrum is 
determined by the sign of the larger (atmospheric) neutrino mass-squared difference. It is 
clear, however, that the sign of the atmospheric mass-squared difference has no physical 
meaning: it is a convention based on the labeling of the neutrino masses and determination 
of the neutrino mass-squared difference (Am^j. = — mf). In both cases of the neutrino 
mass spectrum for the mixing angles the same notations can be used. 
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In the case of the normal spectrum the neutrino masses m2,3 are connected 
with the hghtest mass mi and two neutrino mass-squared differences Amjg 
and Am23 by the following relations 



m2 = Y + Am^2 ) ms = y + Ami2 + Am^s (161) 
In the case of the inverted spectrum we have 



mi = Y m| + I Amfgl, m2 = y m| + | Am^gl + Am^2 (162) 

It is obvious that effective Majorana mass is determined not only by the 
lightest neutrino mass and neutrino mass-squared differences but also by the 
character of the neutrino mass spectrum. 

Usually the following three typical neutrino mass spectra are considerecJ^ 



^^Let us notice that these three neutrino mass spectra correspond to different mecha- 
nisms of neutrino mass generation. Masses of quarks and charged leptons satisfy hierarchy 
of the type (|163|) . Hierarchy of neutrino masses is a typical feature of GUT models (like 
50(10)) in which quarks and leptons are unified. Inverted spectrum and quasi-degenerate 
spectrum require specific symmetries of the neutrino mass matrix. 
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1. Hierarchy of the neutrino masses 



mi <^ m2 <^ ms. (163) 

2. Inverted hierarchy of the neutrino masses 

<^ mi < 1712 (164) 

3. Quasi-degenerate neutrino mass spectrum 



mi ~ 772,2 — ma, ''^i("^3) ^ y ^"^23 ( y I ^"^isl)- (165) 

We will discuss now the possible values of the effective Majorana mass in the 
case of these three neutrino mass spectra. 

I. Hierarchy of the neutrino masses 

In this case we have 



mi ^ y/Amf2, rn2 — yAmfg, rn^ — y Aml^. (166) 

Thus, in the case of neutrino mass hierarchy the neutrino masses m2 and 
are determined by the neutrino mass-squared differences Am^g Am23, 
correspondingly, and the lightest mass is very small. Neglecting the con- 
tribution of mi to the effective Majorana mass and using the standard 
parametrization of the neutrino mixing matrix we find 



\m(si3\ 



sin^ 012 V Amf2 + e^^" sin^ ^13 J Am^g 



(167) 



Here a is a Majorana phase difference. 

The first term in Eq. fll67p is small because of the smallness of Am^2- The 
contribution of the "large" Am^a to |m^/3| is suppressed by the small factor 
sin2^i3. Using the values ffTOTl) and ffTOSl) and the CHOOZ bound ffT09l) . we 
have 



sin^ ^12 y Am^g ~ 2.8 • 10"^ eV, sin^ ^13 -^Am^g < 2.5 ■ IQ-^eV. (168) 

Thus, if the value of the parameter sin^ 6'i3 is close to the CHOOZ bound, 
the first term and the modulus of the second term of f ll67p are approximately 
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equal and at a 7i/2 the terms in the expression (I167p practically cancel 
each other. In this case the Majorana mass will be close to zero. 

Even without this possible cancelation the effective Majorana mass in the 
case of the neutrino mass hierarchy is very small. In fact, from fll67p and 
(11681) we have the following upper bound 



\m,3(3\< (^sin^ Oi2 -\/a^ + sin^ 61^ ^Am^gj < 5.3 ■ 10"^ eV. (169) 

This bound is significantly smaller that the expected sensitivity of the future 
experiments on the search for Oz//3/3-decay (see later). 

II. Inverted hierarchy of the neutrino masses 
For the neutrino masses we have in this case 
ms < ^lAmfgl, mi ~ ^lAm^g], ~ ^lAm^g] (1 + i^^X (170) 



In the expression for the effective Majorana mass |wi/3/3| the lightest mass 
is multiplied by the small parameter sin^6'i3. Neglecting the contribution of 
this term and also neglecting the small term | in f ll70p we find 



Amfgl (1 -sin2 2 0i2 sin^a)!, (171) 



where a is the difference of the Majorana phases of the elements Uf.2 and Uei- 
The phase difference a is the only unknown parameter in the expression for 
|m^/3| in the case of the inverted hierarchy. From fll7ip we find 



cos2^i2 y lAm^gl < \mpp\ < sj\/^mlj\. (172) 

The upper and lower bounds of the inequality (11720 corresponds to the CP- 
invariance in the lepton sector. In fact, the elements of the first row of the 
neutrino mixing matrix can be written in the form Uei = \Uei \ e*"\ In the case 
of the CP- invariance, the elements of the neutrino mixing matrix satisfies 
the condition ( IHTl) . From this condition we have 

e''"' = (173) 



35 



where 77, = ±i is the CP parity of the Majorana neutrino with mass m^. For 
the phase difference a = 02 — «i we have 

e'^" = V2V*i- (174) 

If ri2 = rji we obtain a = 0, vr (the upper bound in the inequahty fll72p ). If 
ri2 = —rji we have a = ±^ (the lower bound in the inequality fll72p ). 

From f ll07p and fllOSp we find the following range of the possible values 
of the effective Majorana mass 

1.8- 10^2 < i^^^l < 49.10^2 eV (175) 

Thus, in the case of the inverted hierarchy of the neutrino masses the lower 
bound of the effective Majorana mass is different from zero. 

The anticipated sensitivities to the effective Majorana mass of the next 
generation of the experiments on the search for the Oz//3/3-decay are in the 
range fll75p (see below). Thus, the future Oz//3/3-decay experiments will probe 
the Majorana nature of neutrinos with definite masses in the case of the 
inverted hierarchy of the neutrino masses. 

III. Quasi-degenerate neutrino mass spectrum 

Neglecting the small contribution of sin^^is, for the effective Majorana 
mass we obtain in the case of the quasi-degenerate neutrino mass spectrum 
the following expression 

|m/3/3| ~ mmin (1 - 8111^2612 siu^tt)^, (176) 

where mmin is the lightest neutrino mass and a is the Majorana phase differ- 
ence. Thus, Imi^pl depends in this case on two unknown parameters: m^-^ 
and a. From fll76p we obtain the following range for the effective Majorana 
mass: 

cos2 6'i2mmin < |m^/3| < (177) 

If Oz//3/3-decay will be observed and the effective Majorana mass turn out to 
be relatively large (|m^^| it would be an evidence that neutrinos 

are Majorana particles and the spectrum of their mass is quasi-degenerate. 
In this case we could conclude that 

\m^f3\ < mmin < 2.8 |m^/3|. (178) 
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An information about the lightest neutrino mass can be obtained from 
experiments on the measurement of the end-point part of the /9-spectrum 
of tritium. From existing data of the Mainz and Troitsk [50] tritium 
experiments it was found the upper bound 

m^nin < 2.2 eV. (179) 

The sensitivity of the future KATRIN experiment [51] is expected to be 

minin ^ 0.2 eV (180) 

We have considered three neutrino mass spectra with special values of the 
lightest neutrino mass mmin- In Fig- 1 the effective Majorana mass for the 
normal and inverted neutrino mass spectra as a function of mmin is presented. 
Uncertainties of the parameters Am^g; ^"^23 ^^"^ ^^^^ ^12 and possible values 
of the Majorana phase difference a are taken into account in Fig.l. 




1 10 100 1000 

Minimum Neutrino Mass (meV) 



Figure 1: The effective Majorana mass for the normal and inverted neutrino 
mass spectra as a function of minimal neutrino mass. 

In conclusion let us notice that if in the KATRIN (or other) experiments 
the neutrino mass will be measured and in the Oz^/3/3-decay experiments, 
sensitive to the effective Majorana mass in the range fll77p . a positive signal 
will not be observed it would be an evidence that neutrinos with definite 
masses are Dirac particles. 
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8 Nuclear matrix elements of Oz//3^-decay 



The effective Majorana mass {rrii^pl is not directly measurable quantity. From 
the measurement of the half-life of the 0i//9/9-decay only the product of the 
effective Majorana mass and the nuclear matrix element can be obtained . 
In order to determine the effective Majorana mass we need to know nuclear 
matrix elements of the Oz//3/3-decay (NME). 

The calculation of NME is a complicated nuclear many body problem. 
Two main approaches are used: Nuclear Shell Model (NSM)[55j and Quasi- 
particle Random Phase Approximation (QRPA)[56| [57]. 

The Nuclear Shell Model is attractive from physical point of view: there 
are many spectroscopic data in favor of shell structure of nuclei (spins and 
parities of nuclei, binding energies of magic nuclei, etc.). It is based on the 
assumption that exist spherically symmetrical averaged nucleon field (usually 
oscillator potential) and one-particle states in this field are used as a basis 
for the description of valence nucleons. An effective interaction between 
nucleons is taken into account in the Hamiltonian. Because of computational 
difficulties, rather limited number of one-particle states can be used in the 
NSM calculations. However, all possible distributions of valence nucleons 
over these states are taken into account. 

The neutrinoless double /3-decay of a nucleus is due to transition of two 
neutrons into two protons with the emission of two electrons. The operator 
of the transition of two neutrons into two protons can be presented in the 
form of the sum of products of a operator of the absorption of two neutrons 
in a state with total momentum J and parity vr and an operator of creation 
of two protons with the same momentum and parity : 



It was found [58] that the dominant contribution to the NME comes from 

the 0^ state of the neutron-neutron pair. Sizable contribution gives also 
2"*" state. It is, however, smaller and has opposite sign. The contributions 
of other states are negligibly small. The dominance of the contribution of 
the 0^ state corresponds to the pairing content of the initial and final wave 
functions. Let us notice that if seniority of the initial and final wave functions 
is equal to zero NME would be maximal. 

Further, it was found [52] that the major contribution to NME comes 
from pairs of neutrons at the distance r < {2 — 3) fm. In order to take 




(181) 
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into account strong repulsion of nucleons at small distances (< 1 fm) addi- 
tional r-dependence (so called short-range correlations) is introduced in the 
expression for the NME. 

This additional r-dependence is parameterized by Jastrow-type function 

m 

f{r) = 1 - e-'^'^'il - hr"^), a = 1.1 fm'^ b = 0.68 fm'l (182) 

Recently it was proposed to take into account the short-range correlations 
by an Unitary Correlation Operator Method (UCOM) |6T]. In this method 
the correlated wave function is obtained by an unitary transformation of 
uncorrelated wave function. 

In the Table [8] we present the NME values of nuclear matrix elements 
of the Oi//3/3-decay M^" which were calculated with Jastrow-like and UCOM 
short-range correlations. 

Table I 

The NSM values of nuclear matrix elements of the Oz//3/3-decay |59] 



Nuclei transition 


MO'^(UCOM) 


M°''(Jastrow) 


48Ca -> ^^Ti 


0.85 


0.64 


™Ge ^ ^6gg 


2.81 


2.30 


82Se ^ 82Kr 


2.64 


2.18 




2.62 


2.10 


i28Te -> i28Te 


2.88 


2.34 


i30Te ^ i30xe 


2.65 


2.12 


i36xe ^ i36Ba 


2.19 


1.76 



Notice that except double-magic nucleus Ca NSM nuclear matrix ele- 
ments of the Oz//3/9-decay for all considered nuclei are practically the same 
(they differ by not more than ~ 20 %). 

There are two groups which are performing QRPA calculation of NME 
at present: Tubingen group ^21 ESI El] and Jyvaskyla group [SSI EEl EH 155] . 
QRPA method allows to include pairing correlations in nuclear wave func- 
tions trough the introduction of quasiparticles (particle- hole pairs). Two pa- 
rameters Qpp and Qph of the model characterize particle-particle and particle- 
hole interactions. The constant Qph is obtained from the fit of the energy 
of the giant Gamov- Teller resonance. The Tubingen group determines the 
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value of the constant gpp from the measured half-hfe of the 2i//3/3-decay of 
the corresponding nucleus. The Jyvaskyla group determines the constant gpp 
from data on the /3-decay of nuclei which are close to the nuclei of the interest 
for the Oz//3/9-decay. They also use the value of the constant gpp, obtained 
from the half-life of the 2i//3/3-decay. 

In QRPA approach the mean nuclear field is described by the Woods- 
Saxon potential. The number of basic one-particle states which can be used 
in the QRPA is much larger than in the NSM. This is an important advantage 
of the QRPA approach. However, only limited excitations can be taken into 
account. 

Like in the NSM case, in the QRPA approach the dominant contribution 
to NME gives O"*" state of neutron pairs. However, in QRPA not only 2"*" 
state but also other states give significant contribution. 

In both approaches main contribution to NME gives neutron pairs at the 
distance smaller than (2 — 3) fm. The short-range correlations, taking into 
account nucleon repulsion at short distances, are introduced in the QRPA 
expression for NME via the Jastrow-type function (11821) and through Unitary 
Correlation Operator Method procedure. Recently [64j the short-range cor- 
relations were calculated directly from different nucleon-nucleon potentials 
by the coupled cluster method (CCM)[69]. 

In the Table II the results of the calculations of the QRPA nuclear matrix 
elements by the Tubingen group are presented. The short-range correlations 
were calculated by CCM method.. For comparison in the Table II the results 
of the calculation of NME with the Jastrow-type short-range correlations are 
also presented. 
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Table II 



The values of QRPA nuclear matrix elements of the Oi//9/3-decay with CCM 
and Jastrow short-range correlations [6l] . 



Nucleus 


M°'' (Jastrow) 




™Ge 


3.33 - 4.68 


4.07- 6.64 


82Se 


2.82 - 4.17 


3.53 - 5.92 


962r 


1.01 - 1.34 


1.43 - 2.12 


looMo 


2.22 - 3.53 


2.94 - 5.56 


looMo 


2.22 - 3.53 


2.94 - 5.56 




1.83 - 2.93 


2.30 - 4.14 


128Tg 


2.46 - 3.77 


3.21 - 5.65 




2.27- 3.38 


2.92 - 5.04 


i36Xe 


1.17 - 2.22 


1.57- 3.24 



The uncertainties of NME in Table II are mainly due to different values of the 
axial constant qa which are used in the calculations. Upper bounds of NME 
correspond to the free nucleon value qa = 1-25 and lower bounds correspond 
to quenched in the nuclear matter value qa = I- 

The results of the calculations of nuclear matrix elements of the Oz//3/3- 
decay performed by the Jyvaskyla group are presented in the Table III. The 
short-range correlations were taken into account by Jastrow and UCOM pro- 
cedure. 
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Table III 



The values of NME calculated in the framework of QRPA by the Jyvaskyla 

group [65] 



\n r* Id! 1 c 

IN UL-lcUo 


ypp 








™Ge 


1.02 


1.00 


5.08 


6.56 




X .uu 




4 03 


5 36 


^^Se 


0.96 


1.00 


3.54 


4.60 




1 nn 




9 7S 


79 


96 Zr 


1.06 


1.00 


3.13 


4.31 




1.11 


1.25 


2.07 


3.12 


looMo 


1.07 


1.00 


3.53 


4.85 




1.00 


1.25 


2.74 


3.93 




0.82(/3) 


1.25 


3.98 


4.93 




0.97 


1.00 


3.68 


4.68 




1.01 


1.25 


3.03 


3.94 


128rpg 


0.86(/3) 


1.25 


4.07 


5.51 




0.89 


1.00 


4.23 


5.84 




0.92 


1.25 


3.38 


4.79 




0.84 


1.00 


4.06 


5.44 




0.90 


1.25 


2.99 


4.22 


i36Xe 


0.74 


1.00 


2.86 


3.72 




0.83 


1.25 


2.05 


2.80 



It is difficult to expect that outcome of the many-body nuclear calculations, 
based on different assumptions, will be the same. However, from the results 
presented in Tables I-III we can conclude the following 

1. The values of the nuclear matrix elements of the Oi//3/3-decay of dif- 
ferent nuclei obtained in the latest QRPA and NSM calculations are 
qualitatively compatible. 

2. NSM nuclear matrix elements of ^^Ge, ^^Se and ^^''Te are by a factor 
(1.5-2) lower than QRPA nuclear matrix elements. 

3. There is no doubts that traditional methods of the calculation of NME 
will be improved and, apparently, new methods will appear. However, 
it will be very important to find a way to test the calculations. 
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If neutrinoless double /3-decay will be discovered and half-live of dif- 
ferent nuclei is measured, from the ratios of measured half-lives in this 
case it will be possible to test different models of the calculation of NME 
[To] . If, for example, half-lives of the Oz//3/3-decay of ^^Ge and ^^°Te will 
be measured, the ratio of half-lives will be practically equal to the in- 
verse ratio of the corresponding phase-space factors in the case of NSM 
nuclear matrix elements and could be significantly different from this 
ratio in the case of QRPA nuclear matrix elements. 

9 Experiments on the search for Oz//3/3-decay 

At present exist data of many experiments on the search for neutrinoless 
double /3-decay. The most stringent lower bound on the half-lives of the 
Oi//3/3-decay of different nuclei was obtained in the Heidelberg-Moscow [1] 
and IGEX [2] experiments, and in the recent CUORICINO p] and NEMO 
[71] experiments. 

In the Heidelberg-Moscow and IGEX experiments two electrons with to- 
tal energy Qpp = 2039 keV which are produced in the 0^ — t- 0^ transition 
^^Ge — >■ ^^Se -|- + were searched for. In the Heidelberg-Moscow experi- 
ment the source (and detector) consist of five crystals of 86 % enriched ^^Ge 
with total mass 10.96 kg. In the IGEX experiment ~ 7 kg of enriched ^^Ge 
was used. Low background level (~ 0.06 counts/ (keV kg y)) and high energy 
resolution ( ~ 3 keV) were reached in the germanium experiments. 

For the half-life of ''^Ge in the Heidelberg-Moscow experiment the follow- 
ing lower bound was obtained [1] 



From this result the following upper bound on the effective Majorana mass 
was inferred: < 0.35 eV. 

In the IGEX experiment it was found [2] 



From this results, assuming different NME, it was found the bound: \mpp\ < 
(0.33 - 1.35) eV 

In the cryogenic experiment CUORICINO [3] the search for the Op^jS- 
decay of ^^°Te was performed. An array of 62 Te02 crystals with a total 



T°;2(^'^Ge > 1.9- lO^V 



(183) 



T°;2(^'^Ge > 1.57- lO^V 



(184) 
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active mass of 40.7 kg was cooled to (8 — 10) mK in a dilution refrigerator. 
Since the heat capacity is proportional to an increase of temperature 
due to tiny release of energy in the Oz//3/3-decay can be recorded by special 
thermometers. 

No evidence for the Oz//3/3-decay of ^'^'^Te was obtained in the CUORICINO 
experiment. For the half-life of ^^^Te a limit 

T°;2('^°Te > 3.0 ■ lO^V- (185) 

was obtained [3]. From this limit using the values of the NME, calculated in 
the latest papers, the following upper bound was inferred |m/3/3| < (0.19 — 
0.68) eV. 

In the NEM03 experiment ^Tlj the cylindrical source was divided in sec- 
tors with enriched i°°Mo (6914 g), ^^Se (932 g), ^^^cd (405 g), ^^^Te (454 
g), i^°Nd (34 g), 96Zr (94 g) and ^^Ca (7g g). For the detecting of the two 
electrons drift cells and plastic scintillator were used. No Oz//3/3-decay was 
observed. In the Table IV the results of the NEM03 experiment are pre- 
sented. 

Table IV 

Lower bounds of the half-lives of the Oi//3/3-decay of different nuclei, 
obtained in the NEM03 experimental]. 



Nucleus 


T0;2 (90%CL) 


\mf3(s\{eV) 


looMo 


^ 5.8 ■ 10^3 y 


^ (0.6-1.3) 


82Se 


^ 2.1 ■ 10^3 y 


^ (1.2 - 2.2) 


96 Zr 


^ 8.6- 10^1 y 


^ (7.4-20.1) 


48Ca 


^ 1.3-1022 y 


^ (29.7 


150Nd 


^ 1.8 - 1022 y 


^ (4.0-6.3) 



Several new experiments on the search for the Oi//3/3-decay are at preparation 
at present. In these new experiments it is planned to reach the sensitivity 
|m^^| ~ a few 10~2 eV , corresponding to the inverted hierarchy of the neu- 
trino mass spectrum. 

In the future GERDA experiment ^72j array of enriched Ge crystals will 
be cooled and shielded by liquid argon (or nitrogen) of very high radiopurity. 
In the phase I of the GERDA experiment 5 detectors from the Heidelberg- 
Moscow experiment (active mass 11.9 kg) and 3 detectors from the IGEX 
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experiment (active mass 6 kg) will be used. The expected background at this 
phase of the experiment will be ~ 10~^ counts/(kg keV y). The expected 
sensitivity will be Ti/2(^^Ge) ~ 3 ■ 10^^ y at 90% CL. Nonobservation of the 
neutrinoless double /3-decay at this phase of the experiment would allow to 
obtain the upper bound \mpf}\ < 0.27 eV (with QRPA NME). 

During the phase II of the GERDA experiment additional 22 kg of the 
enriched Ge will be used (total active mass of the enriched Ge will be about 
40 kg). The expected background 10~^ counts/(kg keV y). The sensitivity 
Ti/2(^^Ge) ~ 1.4-1026 y (at 90% CL) is planned to be reached. This sensitivity 
corresponds to the sensitivity to the effective Majorana mass |m/3^| ~ 0.11 
eV (QRPA NME). 

If goals of Phase I and Phase II will be achieved and the level of the back- 
ground lO"'^ counts/(kg keV y) will be reached it is planned (in collaboration 
with the Majorana Collaboration) to build ~ 1 ton germanium detector with 
the aim to investigate the region of the inverted neutrino mass hierarchy. 

As it is well known, the group of participants of the Heidelberg-Moscow 
experiment claimed that it found an evidence for neutrinoless double /3-decay 
of ^^Ge |73]. For the half-life of the decay the authors obtained the following 
3cr range Ti/2(™Ge) = (1.30 — 3.55) ■ 10^^ y. These values correspond to the 
following range for the effective Majorana mass |m/3^| = (0.24 — 0.58) eV. 
(with NME calculated in [H]). There is no detailed analysis of the systematic 
errors in [73] (see [75]). The only way to confirm or refute the claim is to 
perform more sensitive than the Heidelberg-Moscow experiment (preferably 
^^Ge experiment in order to avoid the NME problem). One of the aim of the 
GERDA experiment is to check the claim made in [73] . 

In the proposed Majorana experiment [76] an array of enriched Ge crystals 
will be installed inside of high purity electroformed copper cryostat. It is 
expected that the background in the Majorana experiment will be a factor 
of 150 lower than in the Heidelberg-Moscow and IGEX experiments. Staged 
approach based of the 60 kg enriched Ge array (60/120/180 kg) is planned. 
The expected sensitivity at the first stage of the experiment (Ti/2('^^Ge) ~ 
5.5 ■ 10^6 y ) will allow to check the claim made in the papers [75] 

In the cryogenic CUORE experiment [77j an array of 19 towers made from 
5x5x5 cm^ Te02 crystals is used as a source (detector). The total number 
of the crystals in the experiment is equal to 988. The total mass of the 
crystals 741 kg of TeOz (204 kg of ^^^Te). In the CUORICINO experiment 
one similar tower of a mass 40.7 kg was used. 

The expected background in the CUORE experiment is 0.01 counts/kg 
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keV y). The expected sensitivity to the half-hfe is Ti/2(^^°Te) ^ 2.5 ■ 10^^ y. 
With the present-day values of NME the following sensitivity to the effective 
Majorana mass will be achieved: Impisl ~ (4.7 — 5.3) ■ 10~^ eV. 

In the future EXO experiment |78] the 0i//3/9-decay of ^^^Xe will be search 
for. Because there is no need to grow crystals and procedure of enrichment is 
relatively simple, Xe is ideal for a large scale (one ton or more) neutrinoless 
double /3-decay experiment. Ion ^'^^Ba'*'^, produced in the decay ^^^Xe — >■ 
^^^Ba^"*" + e~ —, by the capture of an electron can be transferred to the 
ion ^^^Ba"^ which is stable in Xe. The EXO Collaboration plan to identify 
^^^Ba"^ ion by optical pumping with lasers. Single ion can be detected by 
this technique (via photon rate 10^/s). When the program of the ^^^Ba"*" 
tagging will be realized, the background in the experiment on the search for 
the Oi//3/3-decay will be drastically reduced. 

At present the EXO collaboration is constructing 200 kg liquid xenon 
TPC with Xe enriched to 80% in ^^^Xe. No ^^^Ba^ tagging will be done 
at this stage. In this experiment the sensitivity l^/?/?] — 1.5 ■ 10~^ eV is 
anticipated. 

We have discussed experiments on the search for neutrinoless double (3- 
decay which will be done in the coming years. There are several other exper- 
iments which are in R & D stage: Super-NEMO (i5°Nd or ^^Se) [79], MOON 
(io°Mo) iSOJ, SNO++ (i^°Nd)[8T], COBRA {^^^Cd, ^^^Te )[82], CANDLES 
(^«Nd)[M], DCBA (i50Nd)[Sl], CAMEO (ii6Cd)[H5], XMASS (i=^6xe)|HS] and 
others. 



10 Conclusion 

The observation of the neutrino oscillations in experiments with atmospheric, 
solar, reactor and accelerator neutrinos proves that neutrino masses are differ- 
ent from zero and that the states of flavor neutrinos Ue, i/^, utau are mixtures 
of states of neutrinos with different masses. There are two general possi- 
bilities for neutrinos with definite masses: they can be 4-component Dirac 
particles, possessing conserved total lepton number which distinguish neu- 
trinos and antineutrinos or purely neutral 2-component Majorana particles 
with identical neutrinos and antineutrinos. 

It will be extremely important for the further development of the theory 
of the neutrino masses and mixing to answer the fundamental question: are 
neutrinos with definite masses Dirac or Majorana particles? 
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Neutrino masses are many orders of magnitude smaller than masses of 
their family partners, leptons and quarks. This fact tell us that neutrino 
masses and masses of leptons and quarks have different origin. The most 
natural possibility of the explanation of the smallness of the neutrino masses 
gives us the seesaw mechanism of the neutrino mass generation. This beyond 
the Standard Model mechanism connects smallness of neutrino masses with 
the violation of the total lepton number at a large scale and Majorana nature 
of neutrino masses. If it will be established that neutrinos with definite 
masses are Majorana particles it will be strong argument in favor of the 
seesaw origin of neutrino masses. 

Investigation of the neutrinoless double ^-decay of nuclei is the only prac- 
tical way which could allow to proof that neutrinos are Majorana particles. 
This is simply connected with the fact that there are huge number of parent 
nuclei in a source. However, even if neutrinos are Majorana particle prob- 
abilities of the Oi//3/3-decay is extremely small. There are two reasons for 
that 

• The Oz//3/3-decay is the second order in the Fermi constant process. 

• The Oz//3/3-decay is possible due to neutrino helicity-fiip. In the case 
of neutrino mixing this means that the matrix element of the process 
is proportional to effective Majorana mass m/jyj = UlfUi. Small- 
ness of neutrino masses is additional suppression factor in the decay 
probability. 

Experiments on the measurement of the half- lives of such rare process as neu- 
trinoless double ^-decay with severe requirements to background and energy 
resolution are extremely difficult. A big progress was achieved. However, fu- 
ture experiments with about one ton detectors, which will allow to reach the 
region of values of the effective Majorana mass, which is predicted from neu- 
trino oscillation data in the case of the inverted mass hierarchy, is definitely 
a challenge. Taking into account importance of the problem of the nature of 
massive neutrinos, there is no doubts that goals of future experiments will 
be achieved. 

In this review we considered Oi//3/3-decay, driven by the left-handed SM 
weak interaction and Majorana neutrino masses. If total lepton number 
is not conserved and neutrinos with definite masses are Majorana particles 
such mechanism of the Oz^/3^-decay obviously must exist. In the literature 
many other possible mechanisms of the Oi//3/9-decay were considered (for 
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references see, for example, [S]). We shortly discuss here a mechanism due 
to the exchange of a heavy SUSY neutralino. Let us assume that exist a R- 
parity and lepton number violating interaction which induce the transition 
(i — )■ M + e (e is the selectron) . In combination with the standard SUSY 
interaction which induce transition e — > e + x (x is the neutralino) these two 
interactions in the case of the virtual neutralino provide the Oz//3/3 transition 
n + n— T-p + p + e + e. If the constants of the SUSY interactions are of 
the order of the electroweak constant g and if masses of SUSY particles are 
characterized by a scale A in this contribution of these interactions to 

the matric element of the 0z//3/9-decay is proportional to 

MsusY ~ (186) 

This contribution must be compared with the contribution to the matrix 
element of the Oz//3/3-decay of the standard small Majorana neutrino mass 
mechanism 

Mo ~ G% ^ (187) 

Taking into account that Irri/?/?! < 1 eV and < >~ 100 MeV^ we come 
to the conclusion that for A ~ 1 TeV Msusy can be comparable with Mq 
if a hypothetical SUSY interaction which does not conserve i?-parity and 
the lepton number is characterized by the electroweak constant g (for more 
details see [87]). 

I am grateful to Theory Department of TRIUMF for the hospitality and 
to S. Bacca and A. Schwenk for useful discussion of the problem of nuclear 
matrix element of the Oi//3/3-decay. 



A Ettore Majorana 

Great Italian physicist Ettore Majorana was born in Catania (Scicily, Italy) 
on 5.08.1906. His father was an engineer, specialist in telecommunication. 
There were five children in the family!^ 

In 1921 the family moved to Rome. In 1923 E. Majorana finished High 
School and entered the Engineer Faculty of the Rome University. 

Among his fellow-students and friends were E. Segre. and E. Amaldi. In 
1927 Segre and later Amaldi transferred to Physics Faculty and started to 

-"^^For a detailed biography of E. Majorana see E. Amaldi |88] 
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work with E. Fermi who was appointed in 1926 as a Professor of theoretical 
physics at Rome University. 

E. Majorana was famous at Engineer faculty for his extraordinary ability 
of solving difficult mathematical problems. E. Segre convinced E. Majorana 
to meet and to speak with Fermi. At that time Fermi was developing the 
statistical model which is known as Thomas-Fermi model. He explained 
Majorana the model and showed him the table with numerical values of the 
screening potential which he calculated numerically. 

Next morning Majorana returned back to the Institute of Physics with 
his own table of values of the potential. He transformed second order nonlin- 
ear Thomas-Fermi equation into Riccati equation and solved it numerically. 
Majorana and Fermi results coincided. 

A few days later E. Majorana became student of the Physics Faculty. He 
impressed everybody by his lively mind and broad interests. He was very 
critical person. For his criticism he was called in the Fermi group "Great 
Inquisitor" . 

In 1929 Majorana received diploma. His thesis were devoted to the inves- 
tigation of the structure of nuclei and to the theory of the alpha-decay. His 
supervisor was Fermi. 

After doctorate Majorana visited the Institute of Physics for a few hours 
every day. He spend most of his time in library working and studying Dirac, 
Heisenberg, Pauli, Weil and Wigner papers. 

At that time Fermi and his group worked on problems of atomic and 
molecular physics. Majorana wrote six papers on the subject. These papers 
demonstrated profound Majorana's ability of using symmetry properties of 
the states. This allowed him to simphfy the problem and to choose the suit- 
able approximation (which is normal now but was not usual at that time). 
These papers also demonstrated perfect Majorana's knowledge of experimen- 
tal data. 

In 1932 Majorana received teaching diploma ("libero docente"). Commit- 
tee (Fermi, Lo Surdo, Persico) concluded that "the candidate has a complete 
mastery in theoretical physics" . 

In the end of 1931-beginning of 1932 Fermi and his group started to 
concentrate their efforts on nuclear physics. 

After discovery of the neutron by Chadwick (1932) Majorana was one of 
the first who came to an idea that constituents of nuclei are protons and 
neutrons. He started to develop the theory of nuclear forces. Majorana 
proposed the theory of space exchange forces between p and n (Majorana 



49 



potential). 

Fermi was very interested in the idea and tried to convince Majorana to 
publish his results. However, Majorana refused and even did not allow Fermi 
to mention them in his talk at a conference in Paris. E. Fermi managed, 
however, to persuade Majorana to go to Leipzig where W. Heisenberg was 
working and to Copenhagen where N. Bohr, was working. 

E. Majorana was abroad during seven months, starting from January 
1933. Heisenberg, who worked at that time on the theory of nuclear forces, 
discussed with Majorana his paper on nuclear theory. He convinced Majorana 
to publish it. 

After returning from Germany E. Majorana started to come to the Insti- 
tute of Physics at via Panisperna rather rare and after some months did not 
come at all. 

He was at home and became interested in political economy, philosophy, 
construction of ships, medicine. He even wrote a paper on statistical laws 
in physics and social sciences which was discovered and published after his 
disappearance. 

Meanwhile new talented physicists grown up in Italy (Wick, Racah, Gio- 
vanni Gentili Jr. and others). It was time to create a new chair in theoretical 
physics. This chair was created at the University of Palermo and in the be- 
ginning of 1937 a competition for the chair was announced. 

It was a problem to convince Majorana to take part in the competition. 
Finally, Fermi, Amaldi and Segre managed to convince him. 

Majorana had no publications during several years. He sent to "Nuovo 
Cimento" his most important paper " Symmetrical theory of the electron and 
the positron" in which the theory of the Majorana particles was proposed. 

After that the following happened. By the request of Senator Giovanni 
Gentili E. Majorana for his extraordinary abilities without competition was 
appointed as a professor at Napoli University. 

In January 1938 E. Majorana came to Napoh. In Napoli he had rather 
lonely life. He went to the University only when he had lectures (on Quan- 
tum mechanics). After lectures he visited Professor Carrelli with whom he 
became friendly and discuss different problems in physics. He never men- 
tioned what he was doing. He discussed his neutrino theory and Carrelli had 
an impression that Majorana considered this theory as his most important 
contribution to physics. 

On March 23 1938 E. Majorana decided to go to Palermo. On March 25 
Carrelli received a telegram from Majorana from Palermo. He asked Carrelli 
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do not worry about a letter which he would receive. In the letter which came 
soon, Majorana wrote that he found his life useless and decided to commit 
suicide. Carrelli called Fermi and Fermi called to Luciano, Ettore brother. 
Luciano immediately went to Napoli. He understood that on the evening of 
March 25 Ettore took boat to Napoli. He was seen sleeping in his cabin when 
the boat was entering into the Napoli bay. He did not arrive to Napoli. His 
body was never found. 

During several months there was an investigation conducted by family and 
by the police. Vatican tried to find out whether he entered some monastery. 
No traces were found. 

I will finish with two citations: "There are various kind of scientists in 
the world. The second and third-rate ones do their best but do not get very 
far. There are also first-rate people who make very important discoveries 
which are of capital importance for the development of the science. Then 
there are genius like Galillo and Newton. Ettore Majorana was one of these. 
Majorana had greater gifts that anyone else in the world; unfortunately he 
lacked one quality which other men generally have: plain common sense" 
(E.Fermi from Cocconi memories) 

"E. Majorana was very critical to himself and other people. He was 
permanently unhappy with himself. He was a pessimist but had very acute 
sense of humor. He was conditioned by complicated and absolutely nontrivial 
living rules. ...E. Majorana was quite rich and I can not avoid thinking that 
his life might not have finished so tragically should he have been obliged 
to work for a living. For that reason and also because he did not like to 
publish the results of all investigations he had made, Majorana contribution 
to physics is much less than it could be "(B. Pontecorvo [89j). 

In conclusion I will discuss briefly the content of the Majorana paper 
" Symmetrical theory of electron and positron" . 

E. Majorana was not satisfied with the existing at that time theory of 
electrons and positrons in which positrons were considered as holes in the 
Dirac sea of the states of electrons with negative energies. He wanted to 
formulate the symmetrical theory in which there is no notion of states with 
negative energies. 

Let us consider the Dirac equation for a complex field iplx) 

{ij^da - m) ipix) = 0, (188) 
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where m is the mass of the particles-quanta of the field. The conjugated field 

^^(x) = C7/;^(x), (189) 

(C is the matrix of the charge conjugation) obviously satisfies the same equa- 
tion 

{i-f^da - m) V'(x) = 0. (190) 
Let us present the field ip{x) in the form 

== (191) 

where 

Xiix)- , X2{x)- -^^ . (192) 

It is obvious from f llSSp . f ll90p and f ll92p that the fields Xi,2{x) satisfy the 
Dirac equations 

(27°9„-m) Xi,2(x) =0. (193) 
The fields Xi,2{x) satisfy also additional (Majorana) conditions 

xUx) = Xi,2(x). (194) 

Majorana used the representation in which 7" are imaginary matrices (Ma- 
jorana representation). In this representation iIj'^{x) = iIj*{x) and Xi{x) and 
X2{x) are real and imaginary parts of the field iplx). 

Majorana build quantum field theory of the fields Xi,2{x)- First of all it is 
easy to show that there are no electromagnetic currents for the fields Xi,2{x)- 
In fact, taking into account (11930 . we have 

(x) = Ux)rUx) = -xJ{x){rfUxf = -Ux)rXi{x) =0; {t = 1,2) 

(195) 

Therefore, Xi,2{x) are fields of particles with electric charge and magnetic 
moment equal to zero. 

For the operator of the energy and momentum Majorana obtained the 
following expressions 



Pa = J $^P.(a;(p))t alip) d'p (z = 1,2). 



(196) 
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where operators a^(p) and (a* (p))^ satisfy usual anticommutation relations. 

Thus, {al{p)y ( ci*(p)) is the operator of the creation (absorption) of a 
particle with momentum p and helicity r. There are no states with negative 
energies and quanta of the fields Xi,2{x) are neutral particles (which are 
identical to their antiparticles). 

In the case of the complex field ipi^x) = ^^^^ the current ja{x) = 
ipi{x)''y'^ipi{x) is different from zero. After quantization Majorana came to 
symmetrical theory of particles and antiparticles with operators of total mo- 
mentum and total charge given by the following expressions 

= / + ^^(P) Mp)]d'p (197) 

Q = ^ I Y}^l{Pyr{p) - dlip) dr{p)]d?p (198) 

Here c\.{p){cr{p)) is the operator of the creation (absorption) of particle with 
charge e, momentum p and helicity r and d\,{p){dr{p)) is the operator of the 
creation (absorption) of antiparticle with with charge — e, momentum p and 
helicity r. Correspondingly, 

b), = ct(p)|0), \p)-a = dl{pm (199) 

are states of particle with charge e, helicity r and mass m and antiparticle 
with charge — e, helicity r and the same mass m. 

Majorana wrote in the paper [90]: "A generalization of Jordan- Wigner 
quantization method allows not only to give symmetrical form to the electron- 
positron theory but also to construct an essentially new theory for particles 
without electric charge (neutrons and hypothetical neutrinos)". And further 
in the paper: "Although it is perhaps not possible now to ask experiment to 
choose between the new theory and that in which the Dirac equations are 
simply extended to neutral particles, one should keep in mind that the new 
theory is introducing in the unexplored field a smaller number of hypothetical 
entities" . 

Soon after the Majorana paper Racah [91] and Furry [92] proposed the 
methods which could allow to test whether neutrino is Majorana or Dirac 
particle. The so-called Racah chain of reactions 

(A, Z) ^ (A, Z + 1) + e" + z/, p+ (A', Z') {A', Z' + 1)+ e- (200) 
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is allowed in the case of the Majorana neutrino and is forbidden in the case 
of the Dirac neutrino. Of course, in 1937 Racah could not know that even 
in the case of the Majorana neutrino the chain fl200p is strongly suppressed 
due to neutrino helicity. 

In 1938 Furry considered neutrinoless double /3-decay of nuclei 

{A, Z) {A, Z + 2) + e- + (201) 

induced by the Racah chain with virtual neutrinos. 
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